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Q H , A new method for measuring the quantities influenced by bin-to-bin and inside 

bin correlations is presented. This method is essential for large multiplicity and/or 
high density of particles in phase space. The method was applied to the two particle 
correlation functions of e + e~ — >W + W~ events. 



1 Introduction 



The usual calculations of statistical errors for entries in histograms and the use of these 
errors in the fitting procedure can bias the measurements if there are several entries from 
the same event. Traditionally, in the past this problem was ignored. The effect is small for 
low multiplicity events. However, for LEP and especially for the future LHC, RHIC, etc... 
experiments this effect is not small at all. These entries are correlated creating bin-to-bin 
and inside bin correlations. Neglecting these correlations, as we will show below, leads 
to a remarkable underestimation of the errors in the measured quantities and less precise 
estimation of the quantities themselves. 

An effective approach has been proposed in [II. The exact method decribed here 
includes all possible correlation effects. 

The method was applied to the two-particle correlation functions of e + e~ — >W + W~ 
events. 

2 The Method 

In this section a new method is presented for construction of the covariance matrix between 
the bins of the histogram for an unbiased measurement of the fitted quantities. 

The presence of bin-to-bin correlations in two particle distributions is unavoidable. If 
there are N positive tracks in the events, each of them has (N — 1) entries in the two- 
particle density P, contributing to different bins of the histogram. Also, due to the finite 
bin width, the same track can also enter several times in the same bin, which is a source 
of inside bin correlations. 

The method is based on classical statistics. Let us consider the z-th event from the set 
of n events and the two-particle density P which is presented in the histogram h l with 
N p bins. 

The histogram H = Yh=i ^ an d values 

Cjk = HtM - H jM( h k - B h /n){\ + l/n) 

were calculated event by event. Here j and k are the bin numbers for the histograms. We 
do not know the correlations and errors for one event. But we know that the different 
events are uncorrelated. Let us consider bin values of the histogram made for one event as 
an random vector with unknown distribution. We have an uncorrelated ensemble of these 
vectors and hence we can estimate the covariance matrix statistically. It is important to 
note, that this algorithm computes both the "technical" and the physical correlations. 
For all events we have the resulting histogram H for the two-particle density P and 
Vjk — Cjk ' n/(n — 1) covariance matrix for this histogram. Note, that the expresion 
for unbiased estimation of sample covariance was used for the calculation of covariance 
between two values[2]. The diagonal terms Vjj of the covariance matrix are assumed to 
be the estimate for the squares of the error aj of the j-th of histogram H. 
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3 The Test of the Method Using the Simulated 
Events 



The precise measurement of the correlations between particles became important for 
e + e~ ^W + W~ events due to possible large impact of these correlations on the mea- 
sured W mass jsJj. Thus, we applied the method to the correlation functions in fully 
hadronic and semileptonic WW events. 

The correlation function R is used to study the enhanced probability for emission of 
particles. For pairs of particles, it is defined as 

"o(Pl,P2j 

where P(pi,P2) is the two-particle probability density, Pi is the four- momentum of particle 
i, and Po{pi,P2) is a reference two-particle distribution which, ideally, resembles P{pi,p 2 ) 
in all respects, apart from the lack of Bose-Einstein symmetrization. The effect is usually 
described in terms of the variable Q, defined by Q 2 = — (pi — P2) 2 = M 2 {nir) —4ml, where 
M is the invariant mass of the two pions. The correlation function can then be written 
as 

which is frequently parametrised by the function 

R(Q) = 1 + \e- r2Q2 . (3) 

In the above equation, in the hypothesis of a spherically symmetric pion source, the 
parameter r gives the RMS radius of the source and A is the strength of the correlation 
between the pions. 

The method described in section 2 was tested on the JETSET simulation 0]. We 
choose the fully hadronic and semileptonic decay of WW pairs with Bose-Einstein corre- 
lations included using the LUBOEI code. The value of A = 0.85 and r = 0.5 fm was used. 
In case of fully hadronic channel the Bose-Einstein correlations were switched on for all 
pions(full Bose-Einstein correlations). 

The correlation matrices pjk = Vjk/(<7j<7k) for like-sign pairs of WW semileptonic 
channel (refered as (2q) mode) and for fully hadronic channel (referred as (4g) mode), 
computed using 100 000 simulated events for each sample, are shown in Fig. 1 and Fig. 2. 
The correlations between bins for the WW fully hadronic channel (Fig. 2) are larger than 
in the mixed hadronic and leptonic channel (Fig. 1). Thus, the bin-to-bin correlations are 
increased with multiplicity, as expected. Notice that the correlations are nearly indepen- 
dent of Q, which shows why the effective approach of [1 is a good approximation in this 
case. 

For the future analysis the 500 samples of 3000 events each for (4g) channel and 500 
samples of 1500 events each for (2q) channel were simulated. For each of these samples a 
histogram of the correlation function R(Q) was built, using the 25 bins of 100 MeV from 
to 2.5 GeV. The simulated R distributions were normalized to unity in the region Q > 
0.8 GeV/c 2 . We performed a x 2 fit to the R(Q) to the form (3) for each of 500 samles. 

The average values of A and r from these "naive" fits were: 

\ 2q = 0.333 ±0.029, (4) 
r 2q = 0.562 ± 0.033 fm (5) 
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for (2g) events and 



A 4g = 0.416 ±0.013, (6) 
r 4g = 0.561 ± 0.012 fm (7) 

for (4g) events. The statistical errors correspond to the average of the errors for 500 
samples. The "pull" of the fitted values of A and r for the simulated (4g) events are shown 
in Fig. 3. A Gaussian fits gave that the errors in parameres A and r are underestimated 
by a factor 1.20 ± 0.05 and 1.30 ± 0.05 for the (2q) events and by a factor 1.42 ± 0.06 
and 1.53 ± 0.06 for the (4q) events. 

The average values of A and r from the 500 fits using the inverted Vjk matrix (calculated 
for each of 500 samples) were: 

\ 2q = 0.332 ±0.034, (8) 
r 2q = 0.556 ± 0.040 fm (9) 



for (2g) events and 



X Aq = 0.403 ±0.017, (10) 
r Aq = 0.565 ± 0.016 fm (11) 



for (4g) events. The "pull" of the fitted values of A and r for the simulated (4g) events 
are shown in Fig. 4. A Gaussian fits gave o"(a-<a>)/o"a=1-06 ± 0.05 and cr( r _ <r> )/cr r = 1.02 
± 0.05 for (2q) channel, and ct(a-<a>)/o"a=0.96 ± 0.05 and <T( r _ <r> )/cr r =1.08 ± 0.05 for 
(4q) channel. The above values are in a good agreement with unity and thus the errors 
are correctly estimated. 



4 Summary 

A model independent method for measuring the quantities influenced by bin-to-bin and 
inside bin correlations is described. A package, as an addition to HBOOK, was written to 
support this new functionality. The method was tested using the simulated WW events. 
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PYTHIA A = 0.85 




Figure 1: The correlation matrix for like-sign pairs obtained using the simulated WW 
(2q) events. 
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PYTHIA A = 0„85 




Figure 2: The correlation matrix for like-sign pairs obtained using the simulated WW 
(4g) events. 
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PYTHIA full BEC X=0.85 
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Figure 3: (a) Pull function for fitted parameter A using a binned uncorelated least squares 
fit in the simulated samples. A gausian fit is supperimposed as a solid line, b) Same as 
(a) but for the parameter r. 
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Figure 4: (a) Pull function for fitted parameter A in the simulated samples using covari- 
ance matrix technique. A gausian fit is supperimposed as a solid line, b) Same as (a) but 
for the parameter r. 
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